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Abstract. A mixture of two types of super-paramagnetic colloidal particles with long range dipolar inter- 
action is confined by gravity to the flat interface of a hanging water droplet. The particles are observed 
by video microscopy and the dipolar interaction strength is controlled via an external magnetic field. The 
system is a model system to study the glass transition in 2D, and it exhibits partial clustering of the 
small particles [1 . This clustering is strongly dependent on the relative concentration £ of big and small 
particles. However, changing the interaction strength r reveals that the clustering does not depend on the 
interaction strength. The partial clustering scenario is quantified using Minkowski functionals and partial 
structure factors. Evidence that partial clustering prevents global crystallization is discussed. 

PACS. 82.70.Dd Colloids - 64.70.P Glas transition of specific systems 



1 Influence of dimensionality on frustration 

It is well known that the macroscopic behavior of crys- 
talizing systems sensitively depends on dimensionality, as 
demonstrated by two examples: In 2D an intermediate 
phase exists between fluid and crystal, the hexatic phase, 
where the system has no translational order while the ori- 
entational correlation is still long-range |2|3|4| . Such a two 
step melting scenario is not known in 3D. The Ising model 
for ferromagnetics shows a phase transition for 2D and 3D 
but not for ID [5]. For amorphous systems, however, it 
was found in experiments [6], simulations [7J, and theory 
[8] that the glass transition phenomenology is very simi- 
lar in 2D and 3D systems, both in dynamics and structure 

m- 

A subtle difference, the local density optimization in 2D 
and 3D, is the following: in 3D the local density optimized 
structure of four spheres is obviously a tetrahedron. How- 
ever, it is not possible to completely cover space in 3D 
with tetrahedra, because the angle between two planes of 
a tetrahedron is not a submultiple of 360° [10]. The den- 
sity optimized state with long-range order is realized by 
the hexagonal closed packed structure or other variants of 
the fee stacking with packing fraction <j) = tt/VTS « 74%. 
The dynamical arrest in 3D is expected to be enhanced by 
this geometrical frustration, because the system has to re- 
arrange its local density optimized structure to reach long- 
range ordei0. The local geometrical frustration scenario is 



1 It is found in 3D hard sphere systems that this geometrical 
frustration alone is not sufficient to reach a glassy state as 
it cannot sufficiently suppress crystallization [11 12 13 , and 
additionally polydispersity is needed [14]. 



different in 2D. There, the local density optimized struc- 
ture and densest long-range ordered structure are iden- 
tical, namely hexagonal. For the glass transition in 2D 
it is therefore expected that an increase of complexity is 
necessary to reach dynamical arrest without crystalliza- 
tion: in simulations an isotropic one-component 2D sys- 
tem has been observed undergoing dynamical arrest for 
an inter-particle potential that exhibits two length scales, 
a Lennard- Jones- Gauss potential with two minima [T5j . 
Other simulations showed that systems of identical par- 
ticles in 2D can vitrify if the mentioned local geometric 
frustration is created artificially via an anisotropic five- 
fold interaction potential [16J. Alternatively, the necessary 
complexity can be created by polydispersity as found in 
simulations [T7j . 

A bi-disperse system in general is simple enough to crys- 
tallize as e.g. seen from the rich variety of binary crystal 
structures in an oppositely charged 3D coulomb mixture 
[T8] . In the system at hand, partial clustering prevents the 
homogenous distribution of particles and the system crys- 
tallizes locally into that crystal structure which is closest 
in relative concentration |9|19|20j . In this way the system 
effectively lowers its energy with a compromise between 
minimization of particle transport and minimization of 
potential energy. However, that means that the resulting 
structure is not in equilibrium, but in a frustrated glassy 
state. This competition of local stable crystal structures 
prevents the relaxation into an energetically equilibrated 
state, i.e. a large mono-crystal. 

It was found that a binary mixture of magnetic dipoles is 
a good model system of a 2D glass former [6] as the dy- 
namics and structure shows characteristic glassy behavior: 
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when the interaction strength r is increased, the system 
viscosity increases over several orders of magnitude while 
the global structure remains amorphous. 
For all measured interaction strengths the system shows 
no long-range order as probed with bond order correla- 
tion functions [9J. However, on a local scale it reveals non- 
trivial ordering phenomena: partial clustering and local 
crystallinity. 

Partial clustering |l|21j means that the small particles 
tend to form loose clusters while the big particles are 
homogeneously distributed. The heterogeneous concentra- 
tion of small and big particles leads to a variety of local 
crystal structures when the system is supercooled mak- 
ing up the globally amorphous structure. This local crys- 
tallinity obviously plays a key role for the glass transi- 
tion in this 2D colloidal system as it dominates the glassy 
structure [9]. Therefore, the phenomenon of partial clus- 
tering is indirectly responsible for the frustration towards 
the glassy state. In section |4] the details of the clustering 
scenario are explained. The dependence on the param- 
eters accessible in experiment and the relation to local 
crystallinity [9] are discussed in section [5] 

Firstly, the experimental setup is introduced. After a 
brief discussion about origin of partial clustering, a mor- 
phological analysis using Minkowski measures is presented 
to characterize and quantify the effect. Finally, the de- 
pendence on relevant parameters like the average relative 
concentration £ and the interaction strength r will be dis- 
cussed using Euler characteristics and the partial static 
structure factors. 



2 Experimental setup 
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Fig. 1. Super-paramagnetic colloidal particles confined at a 
water/air interface due to gravity. The curvature of the inter- 
face is actively controlled to be completely flat; therefore the 
system is considered to be ideal two dimensional A magnetic 
field H perpendicular to the interface induces a magnetic mo- 
ment m in each bead leading to a repulsive dipolar interaction. 



The detailed experimental setup is explained in [22] . It 
consists of a mixture of two different kinds of spherical and 
super-paramagnetic colloidal particles (species A: diame- 
ter oIa = 4.5 /im, susceptibility \a — 7.4 • 10 -11 Am 2 /T, 
density pA = 1.5 g/cm 3 and species B: oIb = 2.8 /ira, 
Xb = 6.6 • 10 -12 Am 2 /T, pB = 1.3 g/cm 3 ) which are con- 
fined by gravity to a water /air interface. This interface is 
formed by a water drop suspended by surface tension in a 



top sealed cylindrical hole (6 mm diameter, 1 mm depth) 
of a glass plate as sketched in Fig. [TJ A magnetic field H 
is applied perpendicular to the water /air interface induc- 
ing a magnetic moment M = %H in each particle leading 
to a repulsive dipole-dipole pair-interaction. Counterpart 
of the potential energy is thermal energy which gener- 
ates Brownian motion. Thus the dimensionless interaction 
strength r is defined by the ratio of the potential versus 
thermal energy: 
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Here, £ = Nb/(Na + Nb) is the relative concentration of 
small species with Na big and Nb small particles and p 
is the area density of all particles. The average distance 
of neighboring big particles is given by I a = l/y/p A . The 
interaction strength can be externally controlled by means 
of the magnetic field H. r can be interpreted as an inverse 
temperature and controls the behavior of the system. 
The ensemble of particles is visualized with video mi- 
croscopy from below and the signal of a CCD 8-Bit gray- 
scale camera is analyzed on a computer. The field of view 
has a size of 1170 x 870 pm 2 containing typically 3 • 10 3 
particles, whereas the whole sample contains about up to 
10 5 particles. Standard image processing is performed to 
get size, number and positions of the colloids. A computer 
controlled syringe driven by a micro-stage controls the vol- 
ume of the droplet to get a completely flat surface. The 
inclination is controlled actively by micro-stages with a 
resolution of a « 1 /xrad. After several weeks of adjusting 
and equilibration this provides best equilibrium conditions 
for long time stability. Trajectories for all particles in the 
field of view can be recorded over several days providing 
the whole phase space information. 



3 Origin of partial clustering 

The origin of the clustering phenomenon lies in the neg- 
ative nonadditivity of the binary dipolar pair potential 
|1|21| . It is not expected in positive nonadditive mixtures 
like colloid-polymer mixtures or additive mixtures like hard 
spheres. In binary mixtures with additive hard potentials 
in 2D, phase separation was found using Monte Carlo sim- 
ulations [23]. In addition to the negative nonadditivity, 
the relation vbb < vab < vaa of the pair potentials has 
to be fulfilled |1)21| . Why this leads to partial clustering 
can be understood as follows: The negative nonadditiv- 
ity prevents macro-phase separation as the negativity of 
the nonadditivity parameter A = 2a ab — (o~aa + cfbb) 
means that particles are effectively smaller in a mixed 
state (cjij = J °° dr{l — ex.p[—Vij(r)/kBT] are the Barker- 
Henderson effective hard core diameters). Thus, a mixed 
configuration is preferred from this. Additionally, the in- 
equality vab < vaa energetically favors direct neighbor 
connections between different species instead of big parti- 
cles being neighbors. In competition to this, the inequality 
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vbb < vab favors the neighboring of small particles. The 
best compromise is achieved in the partial clustered ar- 
rangement: neighboring small particles that are located in 
the voids of big particles. 

The genuineness of the effect was demonstrated by a com- 
parison of computer simulation, theory and experiment [lj . 
There, statistical evidence for the occurrence of partial 
clustering is provided by the static structure factor. The 
structure factor Saa(&) nas the characteristic shape of 
a one-component fluid. In contrast, the structure factor 
of the small particles Sbb {k) has a dominant prepeak at 
small wave- vectors which is statistical evidence for an in- 
herent length scale much larger than the typical distance 
between two neighboring small particles: it is the length 
scale of the clusters. 

The prepeak provides statistical evidence of an inherent 
length scale in the small particle configurations. However, 
the structure factors do not reveal all details of the phe- 
nomenon. For example, the voids seen in the big particle 
configurations (see Figure [2|) are not reflected directly in 
the features of the partial structure factors. To further 
elucidate the scenario, the effect is now investigated from 
a morphological point of view using Minkowski measures 
as this provides additional quantitative insight. 



4 Morphological analysis 

For low interaction strengths r the system is an equi- 
librated fluid as seen from the purely diffusive behavior 
[6] . Assuming that entropy is maximized it might be intu- 
itively expected that particles form an arbitrarily mixed 
state where small and big particles are evenly distributed. 
However, already the inspection of a single snapshot re- 
veals that this is not the case, and the scenario turns out 
to be more subtle. How the system appears in equilib- 
rium at low interaction strengths r is demonstrated in 
Figure [2] There, a configuration with relative concentra- 
tion £ = 41% at r = 5 is separately plotted for big (left) 
and small (right) particles. Big particles are distributed 
more evenly while the small particles form loose clusters. 
Configurations of big and small particles are related be- 
cause small particle clusters are able to push away the big 
particles and form voids in the big particle configuration. 
This connection becomes obvious in the highlighted re- 
gions where two big clusters of small particles create two 
voids in the big particle configuration. This visual impres- 
sion of the configurational morphology will be quantified 
in detail after a brief introduction of the used tools, the 
Minkowski measures. 

Minkowski functionals provide morphological measures 
for characterization of size, shape, and connectivity of spa- 
tial patterns in d dimensions [24] . These functionals turned 
out to be an appropriate tool to quantify clustering sub- 
structures in astronomy, e.g. from galaxies [25]. They also 
give insight into the morphology of random interfaces in 
microemulsions [26] . 

The scalar Minkowski valuations V applied to patterns P 




Fig. 2. Snapshot of a particle configuration (r = 5, rela- 
tive concentration £ = 41%) where big (left) and small (right) 
particles are displayed separately. The clusters of small parti- 
cles (right) fit in the voids formed by the big particles (left) as 
highlighted by circles for two examples. 

and Q in Euclidian space are defined by three types of 
covariances [25] : 

1. Invariance to motion; 

2. Additivity: V(P\JQ) = V(P)+V(Q) - V(PC]Q); 

3. Continuity: continuous change for slight distortions in 
pattern P. 

It is guaranteed by the theorem of Hadwiger that in d 
dimensions there are exactly d + 1 morphological mea- 
sures V that are linearly independent [24] . For d = 2, 
the three functionals have intuitive correspondence^: The 
surface area, the circumference of the surface area, 
and the Euler characteristic x- 

In two dimensions the Euler characteristic x for a pattern 
P is defined as 

X = S-H (2) 

where S is the number of connected areas and H the num- 
ber of holes. 

Morphological information can also be obtained from par- 
ticle configurations. As configurations only consist of a set 
of coordinates, a cover disc with radius R is placed on each 
coordinate to construct a pattern that can be evaluated. 
The Minkowski measures are then determined for differ- 
ent cover radii i?, leading to a characteristic curve for a 
given configuration, explained in the following (for better 
understanding follow curves in Figure [3]) 
The first Minkowski measure (disc area normalized to to- 
tal area) increases from 0% to 100% for increasing radius 
R with a decreasing slope when discs start to overlap. The 
second Minkowski measure (circumference) increases with 
cover radius i?, reaches a maximum, and then decays to 
zero when all holes are overlapped. The third Minkowski 
measure (Euler characteristic) is very subtle and describes 
the connection of cover discs and the formation of holes. 
It allows the most detailed interpretation of a given con- 
figuration. 

A typical devolution of an Euler characteristic x/iV with 
N particles can be divided in three characteristic parts for 
continuously increasing cover disc radius R: 

2 For d=3 a common set of functionals correspond to: vol- 
ume, area, integral mean curvature, and Euler characteristic. 
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1. For small R the curve is constant at x/iV = 1 (nor- 
malized to the number of particles TV). Discs are not 
touching and thus S is equal to the number of parti- 
cles. No holes are present. As x is normalized with the 
number of points, the Euler characteristic amounts to 
X/N = 1. 

2. With increasing R the curve drops and x/iV can be- 
come negative when cover discs are large enough to 
overlap and holes are formed. Therefore, the number 
of connected areas S decreases and holes are forming 
which further decreases x/N . The minimum is reached 
when discs are connected to a percolating network and 
the maximum number of holes has formed. 

3. For large R the curve starts to raise again because the 
holes are collapsing until the whole plane is covered 
with overlapping discs and x/N ^ for ^ oo. 

This qualitative behavior is typical for configurations 
in 2D. However, the individual morphological informa- 
tion is obtained from specific features in the three regions 
as the onset of the fall and raise, characteristic kinks or 
plateaus, and the slope of the fall and raise. The Minkowski 
measures in dependence on an increasing cover disc pro- 
vide a characteristic morphological 'fingerprint' of config- 
urations and therefore statistical evidence of the clustering 
scenario, complementary to structure factors 



[1121] . The 

statistical noise of the curves is remarkably small com- 
pared to that of structure factors as the whole statistical 
information of a configuration is contained in every data 
point. Thus, even small features in the curves are true ev- 
idence for morphological particularities. 
All three Minkowski measures in 2D, area, circumference, 
and Euler characteristic, are averaged over 100 configu- 
rations for a given temperature and no time-dependence 
was found during this period of about half an hour. The 
curves are plotted separately for both species in Figure 
[3] in dependence on the cover-disc radius. Corresponding 
snapshots are displayed in Figure [4] to illustrate charac- 
teristic radii as indicated. The used sample was strongly 
supercooled {T = 662, £ = 42%), i.e. it was not in equi- 
librium. However, the features found at these high inter- 
action strengths are the same as for low r, where the 
sample is in equilibrium. The high interaction strength is 
used here, as the discussed features become clearer, but it 
is assumed in the following that the conclusions on clus- 
tering are also valid for low interaction strengths r. This 
assumption will be justified when the dependency on in- 
teraction strength is discussed in section [5] 
First Minkowski measure: Area 
The upper graphs of Figure [3] show the area per parti- 
cle in dependence on the cover-radius R/Ia- The covered 
area starts at zero and is increased continuously to 100% 
when the discs completely overlap the area. The solid line 
in both plots indicates the area fraction covered by non- 
touching free discs. The deviation of the first Minkowski 
measure from that line shows how homogeneous a con- 
figuration is. A clear difference is found between parti- 
cle species: The big particles follow this reference line up 
to « 80%. This is close to the maximum possible value 
for hard discs between 84% and 90.7% for random close 
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Fig. 3. The three 'Minkowski 7 measures in 2D ('area 7 , 'cir- 
cumference ', and 'Euler characteristic ') are shown for a sam- 
ple with £ = 42% at r=662, separately for big (left graphs) 
and small particles (right graphs). Configurations were deco- 
rated with discs and the 'Minkowski' measures are calculated 
in dependence on their radius R/Ia- The solid lines correspond 
to the 'Minkowski' measures of a single cover-disc. Marks A- 
E and a-e in the bottom graphs correspond to radii of interest, 
and an example of a representative configuration at these labels 
is displayed in Figure The features of all measures confirm 
the clustering scenario. 



packing and hexagonal close packing respectively |8)27)28| . 
Therefore, the big particle configurations are very homo- 
geneously distributed. However, the small particles curve 
deviates from the free disc reference much earlier at ~ 40% 
indicating that small particles are much less evenly spread, 
i.e. they form clusters. 

Second Minkowski measure: Circumference 

The middle graphs show how the circumference depends 
on the radius. The black reference line corresponds to 
a free expanding circle. The measure of the big parti- 
cles follows this reference line up to R/Ia ~ 0.5 and 
then sharply decreases. Again, this is due to the homo- 
geneous distribution of the big particles. Their cover-discs 
can expand freely up to the maximum possible value of 
Rmax/^A = 1/2 for square order. The measure of the small 
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Fig. 4. The same section (375 x 375 urn 2 ) of a single particle configuration with cover-discs is shown several times in the 
pictographs A-E and a-e with increasing radius. Pictures correspond to certain radii of interest marked by arrows in the corre- 
sponding graphs of the 'Euler characteristic 7 in Figured In the first (second) row, just the big (small) particles of this particular 
configuration are shown. Voids in the big particle configurations are filled with clusters of small particles ( examples are highlighted 
with circles in pictograph 'A 7 and V respectively) . Configurations were obtained from a sample with £ = 42% at F — 662. 



particles deviates much earlier and less steeply as their 
particle density is very heterogeneous, i.e. clustered. 
Third Minkowski measure: Euler Characteristic 

The most detailed information on morphology is obtained 
from the Euler characteristic. It exhibits many features 
related to characteristic structures of the investigated con- 
figurations. For better understanding of the features de- 
scribed in the following, Figure [H shows snapshots of a 
typical section in the used configuration for specific cover- 
radii. The notation A — E for the big particles and a — e 
for the small ones is used in both Figures. In Figure [4^ 
and[4]A two clusters are highlighted. Firstly, we consider 
only the big particles. The Euler characteristic xa per 
particle is 1 as the expanding discs are not touching for 
small radii (mark A). Again, this continues up to a value 
close to R/Ia ~ 0.5. The characteristic deviation of all 
three measures at this same radius states the homogene- 
ity of the big particle distribution. Then, discs touch and 
Xa/Na decreases rapidly because surfaces are connect- 
ing and holes are forming (mark B). The minimum is 
reached at mark C, and the Euler characteristic immedi- 
ately raises because the smallest holes between the trian- 
gular close packed regions collapse as seen in the compar- 
ison of Figures [UG and [4)3. The next holes to collapse are 
those where one small isolated particle is located. There- 
fore, a little kink is visible at mark D since these one- 
particle holes are a little larger than the holes decaying 
at mark C and therefore 'survive' a little longer. When 
they collapse, the Euler characteristic increases rapidly 
to a pronounced plateau. Note, this plateau is the only 
statistical evidence for the voids in the big particle con- 
figuration made up by the clusters: these voids are large 
and thus they 'survive' for a long 'time' resulting in that 
plateau. These voids are not detectable with the other 
Minkowski measures or the static structure factor of the 
big particles |1|21| . Finally, they start to decay at mark 



E, but not suddenly, which shows that they have a distri- 
bution in size. 

The Euler characteristic of the small particles shows the 
complementary picture: Starting with low values of R/Ia 
the characteristic is Xb/Nb = 1 for free disc expansion. 
The first drop at mark b occurs at much lower values 
than for the big particles because small particles in clus- 
ters connect. The subsequent shoulder right next to mark 
b confirms the clustering: small particles inside a cluster 
are now connected, and it needs some further increase of 
disc radius until the clusters themselves start connecting. 
A small second shoulder at mark c originates from the iso- 
lated particles that are not arranged in clusters. They are 
the last particles incorporated until all discs form a per- 
colating network at the minimum at mark d. The increase 
of Xb/Nb shows how the holes are closing. While the in- 
crease in the Euler characteristics of the big particles has 
a plateau at mark E, the small particles have a clear dip at 
mark e. This reveals information on the shape of the clus- 
ters: The voids in the big particle positions are compact 
in shape stopping the increase of the Euler characteristic 
before mark E. In contrast, the small particles arrange in 
chain-like clusters. When the voids between these struc- 
tures close, they decay into several sub-holes causing the 
characteristic to decrease again. In fact, the big particles 
also cause a little dip at their plateau for the voids can 
sometimes also decay into sub- holes. However, this dip is 
much smaller than for the small ones. 
Most features are also visible in the Euler characteristic 
obtained from Brownian dynamics simulation pQ. There, 
the same qualitative behavior is found but the smaller 
features are 'washed out' because the used interaction 
strengths were much lower, as discussed in the following 
section. Further, the variation of the relative concentra- 
tion £ and the subsequent dependence of the features in 
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Fig. 5. Fraction of small particles (relative to all small parti- 
cles Ns, left axis) that are arranged in clusters. This fraction 
is independent of F but dependent on relative concentration £. 
The dashed lines correspond to the average fractions over all 
interaction strengths T. 

Solid squares represent the fraction of small particles (relative 
to all small particles Ns, right axis) that are arranged in square 
symmetry as evaluated in f^f. 



the Euler characteristic will confirm the interpretation of 
the scenario. 



5 Dependence of clusters on interaction 
strength and relative concentration 

In order to demonstrate the connection between clustered 
equilibrated fluid and supercooled local crystalline struc- 
ture, the dependence of partial clustering on the interac- 
tion strength T and on the average relative concentration 
£ is now discussed. 

In Figure [5] the fraction of small particles arranged in clus- 
ters is plotted versus interaction parameter T for two dif- 
ferent relative concentrations £. A small particle is char- 
acterized as 'cluster-particle', when the closest neighbor is 
also a small particle. This simple criteria implies that the 
smallest possible cluster consists of two close small par- 
ticles surrounded in a cage of big ones. In the graph of 
Figure [5] for £ « 43% it is found that a high fraction of 
~ 80% of all small particles is arranged in clusters. Even 
for a lower relative concentration £ « 29%, still « 65% 
are arranged in clusters. Note, that the fraction of small 
particles in both samples is smaller than that of the big 
particles as £ < 0.5. Therefore, every small particle could 
have enough possibilities to arrange far away from the 
next small particle which is obviously not the case. For an 
arbitrary distribution of the small particles over the num- 
ber of possible sites (which is equal to the number of big 
particles) a fraction of 40% is expected for a relative con- 
centration £ = 29% and a fraction of 55% for a relative 
concentration £ = 43% 0. The fact that these expected 





3 A simple simulation is performed where Na sites are ran- 
domly occupied with small particles, and the same analysis to 
determine the number of cluster-particles is applied. 
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Fig. 6. 7 Euler characteristics' of big (left) and small (right) 
particles for different interaction strengths F . Characteristic 
features are visible for the lowest r and become clearer with 
increasing r. The relative concentration was £ ~ 43%. Curves 
are shifted for reasons of clarity. 



values are significantly lower than the actually measured 
ones additionally confirms that small particles effectively 
attract each other and therefore cluster. 

The main result from Figure [5] for the structure of 
this colloidal glass former is that the fraction of cluster- 
particles is independent of the interaction parameter r, 
in contrast to local crystallinity which strongly increases 
upon supercooling: Clusters do not vanish although the lo- 
cal structure is dominated by local crystallinity for strong 
supercooling [9]. This is demonstrated for the example of 
square order in the same Figure [5] (for details see [9]). 
The local relative concentration is frozen in. Small parti- 
cles are not redistributed to match an equilibrium crys- 
tal structure which would reduce the number of cluster- 
particles (e.g. in square order). In fact, the independence 
of the clustering from r shows that the opposite is the 
case: The clusters force the local structure into that crys- 
talline order which matches best with the local relative 
concentration. In this way, local crystallinity is established 
without long-range order [9] as it inherits the clustered dis- 
tribution of the small particles. 
This behavior is confirmed by the graphs of Figure El 
The Euler characteristics for both species are plotted over 
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Fig. 7. Static structure factors and 'Euler characteristics' of 
samples with comparable interaction strength r but different 
relative concentration £. Indicated values for F and £ are valid 
for all graphs. 

Upper graphs: Static structure factors for big (left) and small 
particles (right). With increasing relative concentration of 
small particles the features of SBB(k) gain more contrast, and 
peaks in the Saa (fc) are shifted slightly towards higher k-values. 
Lower graphs: 'Euler characteristics' for big (left) and small 
particles (right). For decreasing £ the minimum in xb/Nb 
is shifted towards larger radii R. The weight of the features 
changes, e.g. the last shoulder in the drop in Xb/Nb (red ar- 
row, right plot). This shoulder results from the incorporation 
of isolated small particles into the network (see Figure^ and 
is only visible for the lowest relative concentration £ = 30% 
(red curve) at this interaction strength r . The blue arrow in 
the bottom left plot marks the onset of the inter-particle con- 
nection of the big particles for the sample with highest relative 
concentration £ = 61%. The onset is shifted to a lower value 
R compared to the other curves. 



a wide range of the interaction strength r, from fluid 
(r = 5) to the strongly supercooled state (F = 662): the 
curves change continuously. The main clustering features 
as discussed in Figure [4] are visible for all values of r, they 
just become sharper with increasing interaction strength. 
The smallest features like the kink at mark D in Figure 
[3] are smeared out for low F but the plateaus and shoul- 
ders characterizing the partial clustering are qualitatively 
independent of F. 

In Figure the dependence of the local structure on 
relative concentration £ is shown using structure factors 
and Euler characteristics. There, samples with compara- 
ble interaction strengths r but different relative concen- 
trations £ are compared. Adding small particles is shifting 



the peaks of the structure factor SUa(&) towards higher 
fc-values. This can be understood by the clustering effect: 
Small particles form clusters and push the big particles 
closer together resulting in a shift of the main peak. This 
shift is small for the used parameters. However, confirma- 
tion of this interpretation is found in [21] where Liquid 
integral equation theory shows the same result unambigu- 
ously. There, parameters were used that are not acces- 
sible in the experiment (different ratios of the magnetic 
moments xb/xa)- The contrast in Sbb (&) is increased 
for higher relative concentrations £ which is also in agree- 
ment with theory [2T] , 

The Euler characteristics for the same samples, shown in 
the lower graphs of Figure [7j confirm this interpretation. 
The drop in the Xa/Na (bottom left) becomes deeper 
when less small particles are present. Then, the distances 
between big particles are less distributed due to fewer clus- 
ters. The increase becomes steeper for the same reason: 
clusters of small particles cause larger voids collapsing at 
higher cover-disc radii. It is remarkable, that the onset of 
the steep drop is earlier for high relative concentrations 
(£ = 61%, blue curve) as indicated by the blue arrow. 
Again, this is caused by the small particle clusters that 
push together the big particles. 

A strong dependence is found in the Euler characteristics 
of the small particles: the onset of the first drop is inde- 
pendent on the relative concentration indicating that the 
local density of particles in clusters is not affected (un- 
like that of the big ones). What significantly changes is 
the depth of the first drop. The more small particles, the 
deeper the drop, because more small particles are arranged 
in clusters. The last shoulder, before the Euler character- 
istic reaches its minimum (marked by red arrow), refers to 
the isolated particles (see section H]). Therefore, at these 
interaction strengths this shoulder is only visible for the 
sample with the lowest relative concentration £ = 30% 
(red curve) which has the most isolated particles (com- 
pare also with Figured]). 

The systematic dependence of Euler characteristics and 
static structure factors on relative concentration £ con- 
firms the interpretation of partial clustering of section |4] 
However, the main result of this section is that the prin- 
ciple occurrence of the effect is independent of the inter- 
action strength: The clustering in equilibrium at low in- 
teraction strengths is therefore responsible for the variety 
of local crystallinity at strong supercooling suppressing 
long-range order [9]. 



6 Conclusions 

On a local scale the system reveals the non-trivial ordering 
phenomenon of partial clustering: the small particles tend 
to form loose clusters while the big particles are homoge- 
neously distributed. The origin of this effect is traced back 
to the negative nonadditivity of the dipolar pair potential. 
The detailed scenario is quantified using Minkowski func- 
tional applied to experimentally obtained configurations. 
Changing the interaction strength r reveals that the prin- 
ciple scenario does not qualitatively depend on the inter- 
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action strength, and, as a consequence, the local relative 
concentration is simply 'frozen' in. However, the strength 
of the effect increases with the relative concentration £. 
The clustering effect together with the missing ability of 
the system to reorganize fast enough into an equilibrated 
state (i.e. extended crystal structure |9|19|20| ) is crucial to 
understand the glass forming behavior of this system: The 
partial clustering leads locally to a heterogeneous relative 
concentration £ which then leads for increasing interaction 
strengths r to local crystallinity [9] without long-range or- 
der. It provides the necessary complexity for glassy frus- 
tration in this 2D system and prevents solidification into 
the energetically preferred crystalline or poly-crystalline 
morphologies |19|20| . 
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